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The data of the Bullet Cluster 1E0657-558 released on November 15, 2006 reveal that the strong 
and weak gravitational lensing convergence K-map has an So offset from the E-map. The observed 
E-map is a direct measurement of the surface mass density of the ICM gas. It accounts for 83% of 
the averaged mass-fraction of the system. This suggests a modified gravity theory at large distances 
different from Newton's inverse-square gravitational law. In this paper, as a cluster scale general- 
ization of Grumiller's modified gravity model, we present a gravity model with a generalized linear 
Rindler potential in Randers-Finslerian spacetime. The galactic limit of the model is qualitatively 
consistent with the MOND and Grumiller's. It yields approximately the flatness of the rotational 
velocity profile at the radial distance of several kpcs and gives the velocity scales for spiral galaxies 
at which the curves become flattened. Plots of convergence k for a galaxy cluster show that the 
peak of the gravitational potential has chances to lie on the outskirts of the baryonic mass center. 
For proper parameter values, we obtain a reasonable isothermal temperature of the main cluster 
and simultaneously ameliorated the shape of the convergence k curve. The temperature given by 
the model is T ~ 13.1 keV, slightly lower than the observational value T = 14.81^2 keV. 

PACS numbers: 95.30.Sf, 02.40.Yy, 98.62.Sb, 98.65.Cw 

I. INTRODUCTION 

It has long been known that the gravitational potentials of some galaxy clusters are too deep to be generated by the 
observed baryonic matter according to Newton's inverse-square law of gravitation [l[ . This violation of Newton's law 
is further confirmed by a great variety of observations. To name a few: the Oort discrepancy in the disk of the Milky 
Way |2|, the velocity dispersions of dwarf Spheroidal galaxies Q, and the flat rotation curves of spiral galaxies 
The most widely adopted way to solve these mysteries is to assume that all our galaxies and clusters are surrounded 
by massive non-luminous dark matter [H[ . Despite its phenomenological success in explaining the flat rotation curves 
of spiral galaxies, the hypothesis has its own deficiencies. No theory predicts these matters, and they behave in such 
ad hoc way like existing as a halo without undergoing gravitational collapse. There are a lot of possible candidates 
for dark matter (such as axions, neutrinos et al.), but none of them are sufficiently satisfactory. Up to now, all of 
them are either undetected or excluded by experiments and observations. 

Because of all these troubles, some models have been built as alternatives of the dark matter hypothesis. Their 
main ideas are to suggest that Newton's dynamics is invalid in the galactic scale. A famous example is the MOND 
6]. It supposes that in the galactic scale, the Newton's dynamics appears as 



(1) 



lim u(x) = 1, lim u(x) = x, 



where a is a constant and the value of which is of order 10 -8 cm/s 2 . Dwarf and low surface brightness galaxies 
provide a good test for the MOND Q. With a simple formula and the one-and-only-one constant parameter clq, 
the MOND yields the observed luminosity-rotation velocity relation, the Tully-Fisher relation [81. By introducing 
several scalar, vector and tensor fields, Bekenstein developed a relativistic version of the MOND [9(. The covariant 
MOND satisfy all four classical tests on Einstein's general relativity in Solar system. Later, Sanders and Noordermeer 
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extended the MOND analysis to a sample of 17 high surface brightness, early-type disc galaxies [10( • All these render 
the MOND to be an unparalleled candidate for the most-probable modified gravity model in the galactic scale. A 
natural question for the MOND is: how about it in the cluster scale? 

Although the MOND successfully reduces the discrepancy between the visible and the Newtonian dynamical mass 
(which is also quantified in terms of mass-to-light ratio) to a factor of 2/3, there still remains a missing mass problem, 
particularly in the cores of clusters of galaxies [ll|. The data release of the Bullet Cluster 1E0657-558 in November 
of 2006 posed a serious challenge for modified gravity theories such as the MOND. 

The Bullet Cluster 1E0657-558 was first spotted by the Chandra X-ray Observatory in 2002 12]. Located at a 
redshift z = 0.296 (Gpc scale), it has exceptionally high X-ray luminosity and is one of the largest and hottest 
luminous galaxy clusters in the sky. A high- resolution map of the ICM gas, i.e. the surface mass density H(x,y), 



was reconstructed by Clowe et al. [13j, [14| in 2006. It exhibits a supersonic shock front in the plane of the merger, 



which is just aligned with our sky. The high-resolution and absolutely calibrated convergence K-map of the sky region 
that surrounds the "bullet" was also reconstructed by Bradac and Clowe et al. in their gravitational lensing surveys 
13" 12 1 . The K-map is evidently offset from the S-map. The peak of the K-map lies on the region of galaxies instead 



of tracing the ICM gas of the main cluster, which makes up about 83% of the total baryonic mass of the merging 
system. 



Clowe et al. [13j, |16|, [17j took it as a direct empirical evidence of the existence of dark matter, while whether the 



MOND could fit the X-ray temperature profiles without dark matter component is still in issue |ll|, Il8l - l2ll | . Using 
their modified gravity (MOG), Brownstein and Moffat partly explained the steepened peaks of the K-map, while 
attributing the rest differences to the MOG's effect of the galaxies j22j ]. 

On the other hand, Grumiller [23| presented an effective model for gravity of a central object at large scales recently. 
To leading order in the large radius expansion, the action of his model leads to an additional "Rindler term" in the 
gravitational potential. This extra term gives rise to a constant acceleration towards or away from the source. The 
scale where the velocity profile flattens is v ~ 300 km/s, in reasonable agreement with the observational data. 

In this paper, inspired by their prominent work, we try to construct a modified gravity model at large distances with 
a modified Rindler potential — a linear potential with an exponential cutoff — without invoking any non-baryonic 
dark matter. This is carried out in a Randers-Finslerian spacetime from the view point of the Zermelo navigation 
problem 24j-|26|. Finslerian geometry is a generalization of Riemannian geometry without quadratic restrictions on 
the line element 27j. It is intriguing to investigate the possible physical implication in such a general geometrical 



background. In fact, precedent work have yielded some interesting results 28 32]. One of them was that we found 
that the weak field approximation of the Einstein's equations in Finslerian spacetime is equivalent to the MOND 28] . 
The work in this paper is a cluster-scale generalization of Grumiller's model and it is ensured that in the galactic 
limit, it agrees with both the Grumiller's model and the observational data. An approximately flattened velocity 
profile predicted by our model makes it qualitatively consistent with the MOND at the distance scale of several kpcs. 
The Newtonian limit and the gravitational deflection of light are particularly investigated and the deflection angle is 
given explicitly. 

We use the isothermal King /3-model to describe the observed S-map of a galaxy main cluster. The convergence 
k is also obtained. We find that the gravitational potential peak does not always lie on the center of the baryonic 
material center. Chances are that it will has a bigger value in the outskirts rather than the center. This is one of 
the distinguishing features of the reconstructed K-map of the Bullet Cluster system. Besides, the gravity provided 
by the baryonic material is somehow "enlarged". It is reasonable to suggest that these results may ameliorate the 
conundrum between the gravity theory and the observations of the Bullet Cluster 1E0657-558. 

The rest of the paper is organized as follows. In Section 2, we construct a modified gravity model in Randers- 
Finslerian spacetime with a modified Rindler potential at cluster scales. In Section 3.1, we give the Poisson's equation 
by which the effective lens potential obeys. In Section 3.2, by making use of the effective lens potential, we obtain the 
convergence R of the Bullet Cluster 1E0657-558. The cross section of the calculated S-map is presented. In Section 
3.3, the isothermal temperature of the main cluster is calculated. In Section 3.4, we investigate the performance of 
the model in the galactic scale. Conclusions and discussions are given in Section 4. 



II. THE ZERMELO NAVIGATION MODEL AND GRAVITY IN RANDERS-FINSLERIAN SPACETIME 



Finslerian geometry is a natural generalization of Riemannian geometry without quadratic restrictions on the metric 
2jJ. It is based on a function F called Finsler structure with the property F(x,Xy) = XF(x,y) for all A > 0, where 
x M stands for position and = dx^ jd/r for velocity {p, = 0, 1,2, ... , n). The Finsler structure F represents the line 
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1 , 
-F' 



(2) 



element of Finslerian space and the metric is given by 33 1 

_d d 

9t * u = dy^dy~» \2 

The gravity in Finslerian spacetime has long been studied since 1970s mmmmss 

Randers space is a special kind of Finslerian geometry with the Finsler structure F defined on the slit tangent 
bundle TM\0 of a manifold M as [H, E3|, 



where 



F(x,y) = a(x,y) + f3(x,y) 



a(x,y) = ^a^(x)yvy" 
p(x,y) = b^x)y» . 



(3) 

(4) 
(5) 



Here, a^ v is a Riemannian metric and is a 1-form. Here and after, if not specified, lower case Greek indices (i.e. 
/i, is, a, ...) run from to 3 and the Latin ones (i.e. fc, ...) run from 1 to 3. 

Zermelo aimed to find minimum-time trajectories in a Riemannian manifold (M, h) under the influence of a wind 
represented by a vector field W [2q ]. Shen 41 1 proved that the minimum time trajectories are exactly the geodesies of 
Randers space, if the wind is time independent. In a previous paper (4^ |. we supposed a Randers-Finslerian structure 
F(x,y) under the influence of a "wind" in the radial direction W = W^dx^ 1 = W r dr, to wit 



A V + W^W V 



A = 1 - h^W^W" 



where = h» u W v and W v = \{h^ v - W»W V ). is the Schwarzschild metric 

2GM 



hijdx l dx :> = I 1 — 



dr 2 + r z d9 2 + r z sin z 9dip z . 



(G) 



(7) 



The explicit form of F(x,y) reads 



Fdr 



\ 



2GM\~ 1 



dr 2 + r 2 d9 2 + r 2 sin 2 8dip 2 + \- 2 W 2 dr 2 - X 1 W r dr 



= \ x- 



1 - 



2GM 



dr 2 + A" 1 (r 2 d9 2 + r 2 sin 2 9dip 2 ) - X 1 W r dr 



(8) 



where the second equation exploits the expression of A in ^ assuming <C 1. The relativistic form of © is given 
l 



as 



Fdr = \ -X 2 1 - 



2GM\ 



dt 2 + A" 2 1 - 



2GM^ 1 



dr 2 + A" 1 [r 2 d9 2 + r 2 sin 2 9dip 2 ) - A _1 W r dr 



(9) 



In Chapter 8 of [? ] , the standard form of the proper time interval of a static isotropic or approximately static isotropic gravitational 
field is given as 

dr 2 = B{r)dt 2 - A(r)dr 2 - r 2 {dd 2 + sin 2 edip 2 ) . 

The field equations for empty space R^y = requires that A(r)B(r) = constant. And the metric tensor must approach the Minkowski 
tensor in spherical coordinates, that is, for r — > oo, A(r) = B(r) = 1. Thus we have 

A(r)B(r) = 1 . 

For the Randers-Finsler metric J3J and JSJ, that is 

x2 / 2GM\ 
aoo = —A ( 1 — I . 
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The geodesic equations of Randers-Finslerian space is given as 2 

d 2 x^ 



dr 2 

where 



+ 



(ru a + i tl k\ a )y v y a = o, (io) 



e = — b i = — - V 6 (11) 

and is the Christoffel symbols of the Riemannian metric 3 . Given the Randers-Finslarian structure F(x, y) in 
© , the non- vanishing components of the geodesic equations (fTU|) give rise to the relation between the radial distance 
r and the angle tp of the orbits of free particles, to wit 

1 dr\ 2 ( E \ 2 A / 2GM\ 



r 2 dip J \ JX ) r 2 \ r ) ' ^ ^ 



where E and J are the constants of motion [42f . Introducing a new quantity 

GM 



(13) 



the equation (JT2J) can be rewritten in terms of u as 



du\ 2 ( EGM\ 2 



dp) \ JX J 
Suppose the Finslerian parameter A is given as 



- Xu 2 (l - 2u) . (14) 



GM ( r , 

A = 1 1 + — e . (15) 



r„ \ r, 

1-1 



This will result in an effective Newtonian potential' 

GM GM ( r 



>M = 



r, 



1 + — )e~. (16) 



The first term in (|16l) is the usual Newtonian potential and the last linear term with an exponential cutoff is novel. 
The effective acceleration au has two terms, also 

GM GM r _j_ 
a M = 5 5- • — e r ' . (17) 

At sufficiently large distances, the second term may become dominant and provides a linear acceleration towards the 
source. 

Integrating the equation (|14j) . one obtains the deflection angle of light or in a modified Rindler potential in 
Randers-Finslerian spacetime, to wit 

orW = — /( r ; r s . r e), (18) 

r 

where 

f(r;r s ,r e ) = l + — -j^=± ^ J . r2 V ^ dr' . (19) 



This integration can be computed numerically. The model parameters r s and the cutoff scale r e depend on the specific 
gravitational system and are to be determined by observations. For r ^> r e , <Pm and or — > iG r AI . This is 

what we expect in general relativity and the Newtonian limit. The plot of f(r) for different values of r s and r e are 
presented in Figure [TJ 



We just consider the case that the in (J3j is a closed 1-form, i.e. d/3 = 0. 



3 z,^ = ifitis I Ba S v _ da va . da a 

7 va — 2 " ^ dx a dx s ^ dx< 

4 See the Appendix for details. 
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FIG. 1: Plot for the dimensionless Finslerian factor /(r) in the equation (|19p vs. the radial distance r in unit of kpc. Different 
parameter values are plotted for comparison. The black solid line denotes the best-fit result of our model to the K-map of the 
Bullet Cluster 1E0657-558, as shown in Figure H 



III. COMPARING WITH THE OBSERVATIONS 



A. Deflection Angle hr and Convergence k 



In this section, we use the modified gravity model to calculate the convergence k of the Bullet Cluster 1E0657- 
558. First, we get the effective lens potential in Randers-Finslerian spacetime. Hereafter, we use natural units in 
calculations, i.e. setting the speed of light c = 1. 

Einstein's general relativity predicts that a light ray passing by a spherical body of mass M at a minimum distance 
£ is deflected by the angle 



4GM c _ nr —^ 
a= — - — , s = v af +2/ • 



(20) 



The mass of the lens M can be given as 



M(0 = 2tt / £(0*'de' 
Jo 



(21) 



where £(£') is the surface mass density distribution. It results from projecting the volume mass distribution of the 
"lens" p(r) onto the lens plane (i.e. the (x, y)-plane) which is orthogonal to the line-of-sight direction (i.e. the 
z-direction) of the observer, to wit 



p(r)dz 



(22) 



where z = yjr 2 — x 2 — y 2 = y/r 2 — £ 2 and z out = \J r 2 ut — £ 2 . r out denotes the outer radial extent of the galaxy 
cluster, which is defined as when p drops to p(r out ) ~ 10~ 28 g/cm 3 . 
The "Einstein angle" ((2U)) can be rewritten in a vector form as [43| 



«=4G/ d^mi- u <;l> 



(23) 
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where 

d 2 ? = r dip [ d£f (24) 
Jo Jo 

is the surface element of the lens plane. 

With 6 = -g-, one can easily check that ([23)) satisfies 5 (see Section 4.1 in [441 ]) 

& = Sts^ 9 ^ = V « ^ ' (25) 

where 

^(0 = -±r f E(f ) HI*- f'| d 2 f' , S c = . • (26) 

ttE c J R 2 AnGD h D LS 

S c is the critical surface density of the lens. D$ is the angular distance between the observer and the source galaxy, 
i.e. the background. Z?l is the angular distance between the observer and the lens, i.e. the Bullet Cluster 1E0657-558, 
and Dls denotes the angular distance between the lens and the source galaxy. The lens potential i[>(£) obeys the 
two-dimensional Poisson's equation 6 

A^V 2 ^4, A^l^-(^) , (27) 



In astronomy and astrophysics, the quantity ^j- in the equation (1271) is defined as the convergence K, which is also 
called the scaled surface mass density, i.e. 



k = — . (28) 



Suppose a lens potential 



<Hl) = r s , r e ) = -L f E(?)/(£ r s ,r e ) ln|f- f | d 2 f , (29) 



where 



f(£;r s ,r e )= / f(r;r s ,r e )dz (30) 

J—Zout 

and f(r;r s ,r e ) is given by (fl~9|) . For the inner of the lens system, we have £ = Thus, the potential (|39[) can be 
rewritten as 

m = / S(f)/(l > ;r.,r e )InH'-f|d a |» (31) 

= 4" / ln l£ - ?l ^ ■ (32) 
Given the potential (|32[) and using the equation (|25[) . one can reproduce the deflection angle «r in the model, i.e. 

an(0 = ^-f&r s ,r ). (33) 



5 In the two-dimensional polar coordinates, V.-E ^ = et|r 

6 In general, the Laplacian A in polar coordinates is given as 



19 f 8 \ 1 <>- 



For a (/j-independent V(£)> one nas = 0j an d 
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FIG. 2: The E-map from X-ray imaging observations of the Bullet Cluster 1E0657-558, November 15, 2006 data release [I3.[l7|. 
(a) The entire E-map is presented in the equatorial coordinate system J2000. DEC in the j/-axis is short for "Declination" 
and the RA in the i-axis is short for "Right Ascension" . The bright Shockwave region at the right half of the map is the ICM 
gas of the subcluster. The main cluster gas locates at the brightly glowing region to the left of the subcluster gas. The released 
E-map has 185 x 185 pixels and a resolution of 8.5 kpc/pixel. (b) A subset of the E-map on a straight-line connecting the peak 
of the main cluster to that of the subcluster. The peak of the main cluster is taken to be the referential center of the system, 
i.e. £ = . The peak of the subcluster is located at £ ~ 398 kpc. 



B. The E- and K-Map of Bullet Cluster 1E0657-558 



The E-map reconstructed from X-ray imaging observations of the Bullet Cluster 1E0657-558 is shown in Figure 
There are two distinct glowing peaks in Figure [2^ - the left one of the main cluster and the right one of the 
subcluster. A subset of the E-map on a straight-line connecting the peak of the main cluster to that of the subcluster 
is shown in Figure [2Jd. 

For the Bullet Cluster system, the volume mass distribution of the ICM gas of the main cluster p(r) is phenomeno- 
logically described by the King /3-model 45] [46] 47] 



p(r) = Po 



, r 
1+ — 



-, -3/3/2 



r = \/ x 2 + y 2 + z 2 = Vs+^ 



where the parameters po, r c and j3 are determined to be [22j 

po = 3.34 x 10 5 M /kpc 3 , 

P = 0.803 ±0.013 , 

r c = 278.0 ±6.8 kpc. 

Alp, denotes the mass of the sun. 



The outer radial extent of the Bullet Cluster system is given as 



r ou t 



Po 



-2/3/3 



1 



1/2 



2620 kpc. 



10~ 28 g/cm 3 

The radius of the main cluster is ~ 1000 kpc, thus we have £ = The potential (1321) now becomes 



ttE 



(34) 



(35) 
(36) 
(37) 



(38) 



(39) 



c JR 
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where the effective surface mass density is defined as 



p(r)f(r;r s ,r e )dz , z a 



e = v/2620 2 - e kpc . 



(40) 



Making use of (fTTj)) . (j2"5)l . (|34l) and (|4U1) . one finally obtains the convergence K-map of the Bullet Cluster system 



_ s(0 



Po 





1 + 








l—Zo-M 







-3/3/2 



f(r;r s ,r e )dz 



where 



f(r;r s ,r e ) = 1 + 



1 

2r, 



(2 + 7^) (l+fe) ej-3(l+j 



2(1-^) 



-cZr 



(41) 



(42) 



and 



• the parameters p 0j r c and /? are given in ([55)1 to (1571) . 

• z = \j r 1 — x 2 — y 2 = \/ r 2 — £ 2 and z out is given by (1401) , 

• for the Bullet Cluster 1E0657-558, one has Dl ^ ls ~ 540 kpc. So S c in flU]) takes a value of 

Ds 



4irGD L D LS 



3.1 x 10 y M /kpc z , 



(43) 



• r s and r e are model parameters to be determined by fitting (|41[) to the K-map reconstructed from the gravitational 
lensing survey. 

The K-map obtained from the strong and weak gravitational lensing survey of the Bullet Cluster 1E0657-558 is 
presented in Figure One can see that the two distinct glowing regions in Figure |3k - the left one of the main 
cluster and the right one of the subcluster - somewhat depart from those shown in Figure A subset of the K-map 
on a straight-line connecting the peak of the main cluster to that of the subcluster is also shown in Figure ^p. 

A section of the K-map (|41l) . which crossing the two peaks is plotted in Figure [4] Sections of k- and S-map 
reconstructed from the strong and weak gravitational lensing survey and the X-ray imaging observations are also 
shown in the same graph for comparison. Our approach follows a sequence of approximations: 

• Take the main cluster thermal profile to be isothermal. 

• Neglect the subcluster for zeroth order approximation. 

• Perform the fit using a section of the K-map on a straight-line connecting the peak of the main cluster to that 
of the subcluster and then extrapolating it to the entire map. 

• Take the S-peak of the main cluster as the center of the gravitational system, and project the section of the 
K-map onto that of the S-map to make the two overlay for comparison. 



C.The Isothermal Temperature Profile 

The collisionless Boltzmann equation of a spherical system in hydrostatic equilibrium reads 

t PW<7 + [al - aj^) = -p (r — ^ , (44) 

ar r ^ ar 

where $(r) is the gravitational potential of the system and a r and ag^ are respectively the mass- weighted velocity 
dispersions in the radial and (9, cj>) directions. Given an isotropic gas sphere distribution p(r) with a temperature 
profile T(r), one has 

2 2 k B T(r) 

° r = a ^ = "7^7 ' (45) 
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FIG. 3: The K-map reconstructed from the strong and weak gravitational lensing survey of the Bullet Cluster 1E0657-558, 
November 15, 2006 data release [H, Q3]- ( a ) The entire K-map is presented in the equatorial coordinate system J2000. DEC in 
the y-axis is short for "Declination" and the RA in the 2>axis is short for "Right Ascension" . The bright blurred region at the 
left half of the map illuminates the convergence of the main cluster, while the smaller glowing one to the left corresponds to 
that of the subcluster. The released K-map has 110 x 110 pixels and a resolution of 15.4 kpc/pixel. (b) A section of the K-map 
on a straight-line connecting the peak of the main cluster to that of the subcluster. The peak of the main cluster is located at 
£ ~ —180 kpc and that of the subcluster is located at £ ~ 522 kpc. The £ = point is chosen to be the same with that of the 
E-map in Figure [3}5. 



where k B is Boltzmann's constant, pa — 0.609 is the mean atomic weight and m p is the proton mass. The equation 
flU} becomes 



d_ f k B T(r) 
dr \ fiArrip 



P(r] 



-p(r) 



d$(r) 
dr 



(46) 



For the main cluster of the Bullet Cluster system, the ICM gas distribution p(r) is fit by an isotropic and isothermal 
King /3-model with the temperature T(r) = T. Solving the equation f4"6")l for the gravitational acceleration, one 
obtains 



a(r) 



d$(r) 
dr 



dH P (r)) 



dln(r) 



k B T 
PAm p r 
3/3fc B T 
P Atrip \r 2 + r ( 2 



(47) 



Replacing a(r) in (|4"7|) with the effective acceleration au in (TT71) . to wit 



a(r) = a M {r) 



GM 



rir s 



(48) 



we obtain the relation between the dynamical mass M as a function of the radial position r and the temperature T, 
to wit 



M(r) 



3(3k B T 
PArripG \r 2 + r\ 



(49) 
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FIG. 4: (a) Cross sections of the model-predicted K-map and the E-, Ac-map reconstructed from the November 15, 2006 data 
release [l3l . [l7| . The lines denotes the sections of the K-map ()4ip predicted by the Randers-Finslerian model with a modified 
Rindler potential fl6)l. The model parameters (r s ,r e ) are taken to be (23, 187), (32, 204), (45, 222), and (72,272) in unit of 
kpc, of which the isothermal temperature given by the King /3-model is T ~ 13.1, 12.4, 11.3 and 10.6 keV, respectively. The 
sections of the E- and jt-map obtained by observations are respectively represented by small black dots and circles as in Figure 
[5p and|3p. (b) Plot for the dimensionless Finslerian factor f(r) in (|19|l vs. the radial distance r in unit of kpc. Corresponding 
to the Figure (a), the parameter values of (r s ,r e ) are taken to be (23, 187), (32,204), (45,222), and (72,272) in unit of kpc, 
respectively. 
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On the other hand, the mass profile of the main cluster is given by the King /3-model as 

M(r) = 4tt f p(r')r' 2 dr' 
Jo 

-3/3/2 

r' 2 dr' . (50) 

By equating the two M(r) in (|50|) and (|49|) . say at the radial distance r — 1000 kpc, one can determine the isothermal 
temperature T. The isothermal temperature of the main cluster given by our model is shown in Figure SJa) . 

The detection in X-ray by the Einstein IPC, ROSAT and ASCA observations constrained the temperature of the 
main cluster to be T = 17.4±2.5 keV (with 12.3% error) [H] and T = 14.5^ keV (with 6.5% error) It was later 
reported by Markevitch [HI that T = 14.8±^ keV (with 4.5% error). The center value of the main cluster isothermal 
temperature given by our model is T ~ 13.1 keV, albeit slightly lower, is in agreement with the experimental values. 




D. The Galactic Regime 

The particular function form (1151) of the parameter A is inspired by Grumiller's work 23]. The effective potential 
in his paper was given as 

GM , . 

4>m = + Dr, (51) 



where D is constant and the linear term Dr is called the Rindler acceleration term. The effective potential ([T6|) , which 
corresponds to the Finslerian parameter A in equation (14), is merely a generalization of this one. A more general 
form of equation (|5ip can be written as 



GM . . 

hi = + f{r) , (52) 



where f(r) is a function of the distance scale r. In this paper, f(r) takes a form as 



GM 



r 



r, 



f(r) = 1 + - e-~. (53) 



It is a linear potential like the Grumiller's, but an exponential cutoff is added to avoid possible divergence at large 
distances. Grumiller's potential would not suffer from such a disaster because he only discussed the galactic physics. 
While we try to extrapolate the potential (I5T1) to the cluster scale, we do need to consider this problem. 

One should note that the effective potential (l53l) or (fT6f is not "ad hoc" as it seems. It can recover some features of 
the galactic rotation curves predicted by Grumiller's model, which was considered to be a good phenomenological fit 



to the observational data [23[. From v ~ y/wr and (fT7j) . we obtain a new formula for the velocity profile of a galaxy: 



' GM GM 



r 



r, 



2 



v(r)~\ + - e~. (54) 



To describe galaxies, we assume that the total mass M ~ lO n M0 (instead of M ~ 10 14 M Q for the Bullet Cluster 
system). The plot of profile (|54[) for r s ~ 1 kpc and the cutoff scale r e ~ 80 kpc is shown in Figure [SJ From the figure, 
we can see that our model in galactic limit yields an approximately flattened rotation curve of spiral galaxy. It is 
qualitatively consistent with the MOND and Grumiller's model. The velocity scale where the rotation curve flattens 
is ~ 240 km/s, which is in reasonable agreement with Grumiller's prediction and the observational data. A possible 
divergence of the velocity (|54|) at large radial distances is reconciled by the exponential factor to yield a physical 
result. 

In fact, the MOND's results can be incorporated into a Finslerian geometrical background. In our previous work [28j |. 
we have obtained a modified gravity model as the weak field approximation of the Einstein's equations in Finslerian 
spacetime. A specific Finslerian structure would land our gravity model back at the MOND. 
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FIG. 5: Rotation curves of a spiral galaxy v(r) vs. r in unit m/s vs. m (1 kpc ~ 3 x 10 19 m). The dashed line denotes the 
velocity profile predicted by Grumiller (which is qualitatively compatible with the MOND at the distance scale of several kpcs 
and considered a good phenomenological fit to the observational data [23| ). The solid line denotes the results of our model for 
the same total galactic mass. The dotted line which sinks into the bottom is given by Newton's theory, which fails to account 
for the observations. 



IV. CONCLUSIONS AND DISCUSSIONS 



As a cluster scale generalization of Grumiller's gravity model, we presented a gravity model with a generalized 
linear potential. The modified gravity model can be identified with the MOND. The galactic limit of the model 
shared some qualitative features of Gumillcr's result and the MOND. It yielded approximately the flatness of the 
rotational velocity profile at the radial distance of several kpcs. It also gave observation-compatible velocity scales for 
spiral galaxies at which the curves become flattened. 

We also studied the gravitational deflection of light in such a framework and the deflection angle was obtained. The 
modified convergence re formula of a galaxy cluster showed that the peak of the gravitational potential has chances 
to lie on the outskirts of the baryonic mass center. For the Bullet Cluster 1E0657-558 system, the later refers to the 
center of the ICM gas profile of the main cluster. We gave a reasonable isothermal temperature of the main cluster 
and simultaneously ameliorated the shape of the convergence re curve. Proper system-dependent parameter values 
will render the model capable of describing other normal clusters. A full Bayesian analysis is necessary for drawing a 
more confirmatively quantitative conclusion. 

Another feature of the reconstructed K-map of the Bullet Cluster is the apparent spherical symmetry breaking of 
the system. Both the re-peaks of the subcluster and the main one are shifted outwards from the center of the system 
for ~ 200 kpc. A possible way to explain this phenomenon is to postulate it as the nonlinear effect of high-order 
terms of the final modified gravity model for merging clusters. In certain conditions, this high-order effect becomes 
dominant to generate the right half "hill" and the "saddle" between the two peaks in the re-map Fig{3j Specific 
dynamical mechanism remains to be explored. 
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Appendix 

By combining the non- vanishing components of the geodesic equations (jlOp . one obtains the relation between the 
radial distant r and the time t [42j |. 

AE 2 fdrX 2 . J 2 A E^__ c 



B 2 \dt J r 2 B 

where A(r) = A~ 2 (l - and B(r) = A 2 (l — 2^M). £ is an integration constant (see [H| for details). For 

photons, the constant C = 0. The above equation can be rewritten as 

A ,fdr\ 2 J 2 A 1 

\~dt) + ^E~ 2 -B= - < B5 > 



, dt 

order of these quantities (remembering that the leading order terms of A and B are 1), the equation (1551 becomes 



In the Newtonian limit and the weak-field approximation, the quantities ^y, (^l) ,E 2 — I, are small. To first 



dr\ 2 J 2 A 1 



, dt J r 2 B 
Redefining A in ((TSJ) as A = 1 - ^(1 + j-)e~~^ = 1 + 0a, one has 

1 -A"' 1 



, 0. (56) 



r 

1 1 



(1 + A ) 2 1-^ 

(1-20 A ) fl + — J 



1 + 2^M-2 0A 
r 

= -1 + 20 M , (57) 
where 0m = {4>n + 0a) and 4>n = —^r- is the Newtonian potential. Substituting (|57p back into ([56]) . one obtains 

1 /dr\ 2 J 2 A , 1 

2UJ + ^ + ^ = 2' (58) 



where the effective Newtonian potential 0m is given as 

GM GM / r 



r, 



(59) 




r [kpc] 



